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1, INTRODUCTION 

Massive Multiple-input Multiple-output (MIMO) is considered as a key technology in future 5G 
systems to improve the spectrum efficiency. SG made use of large number of adaptive (smart antennas) at the 
base station (BS) [1]. The basic idea of beamforming is to focus energy towards a receiver. Traditionally 
access points were equipped with omni-directional antennas, named so because they send energy in all 
directions. Since they radiate waves in every direction, receiver antennas need not track every client but it 
comes with its own disadvantages. Alternatively, energy can be focused toward the particular receiver which 
is the principle of beamforming [2]. 

Beam forming is the process of combining the weighted of received signals on an array of sensors 
to improve the directionality. Adaptive beamforming is the ability of the beamformer to receive the signal 
only from the desired direction and to reject all other signals from undesired directions. The weight vector for 
the adaptive beamformer continuously changes based on some adaptive algorithm [3, 4]. 

The beamforming technique is used in smart antennas for transmitting and receiving signals in 
massive MIMO systems. Beamforming tries to produce patterns from the antenna’s beam directed towards 
the target mobile terminal, while neglecting interference signals [5, 6]. Adaptive beamforming is more 
suitable for massive MIMO systems than fixed beamforming because of its ability to eliminate interference 
and reduce power consumption [5]. 

There are two types of adaptive algorithm are Blind and Non blind algorithm. Non blind adaptive 
algorithms require the statistical knowledge of the transmitted signal in order to converge to a weight 
solution, such as LCMV (Linearly Constrained Minimum Variance), and MVDR (Minimum Variance 
Distortionless Response) [7]. LCMV uses power constraints to selects optimum weight vectors to reduce the 
filter’s response. LCMV are inappropriate for application to massive MIMO systems because of its low 
convergence rate. 
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MVDR Beamforming algorithm selects the antenna elements weight while keeping a constant gain 
level of beam steered toward the desired direction. The properties of LCMV and MVDR indicate that MVDR 
is more suitable for massive MIMO systems than LCMV. MVDR is more convenient algorithms for massive 
MIMO systems [8]. 

The Non blind based adaptive algorithm are based on minimization of mean square error between 
received signal and reference signal, such as LMS (Least mean square) for smart antenna system. LMS uses a 
gradient based method of steepest decent and requiring reference signal and no knowledge channel [7, 9]. 
Beamforming could be classified into analogue beamforming, digital beamforming, and _ hybrid 
analogue/digital beamforming. Hybrid analogue/digital beamforming has been developed for massive MIMO 
systems to obtain the advantages of analogue and digital beamforming [10]. 

5G requires a high frequency band signals for achieving an extremely high data rate. To make it 
possible mm-wave beamforming which is considered a wideband beamforming can be implemented [11]. 
Hybrid beamforming, which has emerged a low complexity and cost effective solution for millimeter wave 
communications, are the most promising approach in massive MIMO systems [11]. 

The optimal algorithms that satisfy the requirements of next-generation wireless communication 
systems can be one of the adaptive beamforming, such as MVDR, or a combination of two algorithms, to 
estimate accurately the DOA in terms of the azimuth and elevation angles (2D-DOA) [5]. Many researchers 
proposed beamforming methods that work with massive MIMO. 

ALI et al. [5] classified the optimal beamforming techniques that are used in wireless 
communication and suggested an optimal beamforming technique for massive MIMO systems that satisfy the 
next-generation requirements. 

Andreas F. et al. [10] made a survey of the various beamforming methods proposed in the literature 
and indicates that hybrid beamforming are the most promising approach for reducing the hardware cost and 
training overhead in massive MIMO systems. 

S. Shahab et al. [13] evaluated the performance of MVDR by mean of beam pattern of two 
incoming signals and Signal to Interference plus Noise Ratio (SINR). The results showed that MVDR 
performance improves with the increasing number of elements. Also, the beam pattern depends on two 
factors: the number of elements and the spacing between array elements. The best inter element spacing 
obtained is 0.5 to avoid grating lobes and mutual coupling effects. 

Noh et al. [14] demonstrated the use of hybrid beamforming, which are a mixture of digital and 
analogue beamforming, in massive MIMO systems. 

In this paper a new hybrid method solution for massive MIMO was proposed. This hybrid method 
combines MVDR algorithm for its well-known null steering with LMS algorithm that don’t require adjusting 
the weights of the array. 

In this paper a beamforming algorithm was proposed. This method combines the adaptive blind 
beamforming algorithm, MVDR, with the non-blind beamforming LMS. The next two sections contain 
description of these two algorithms. 


2. RESEARCH MEHTODS 
2.1. Minimum Variance Distortionless Response (MVDR) Beamforming 

The basic idea of the Minimum Variance Distortionless Response (MVDR) algorithm or Capon 
beamformer is to estimate the beamforming coefficients in an adaptive way by minimizing the output power 
of the beamformer under a single linear constraint on the response of the array towards the desired signal 
[13, 15]. The MVDR beam former does not require the knowledge of the directions of the interferences for 
weight vector calculation. It requires only the direction of the desired signal. MVDR weight vector is given 
by: [16] 


-1 
R 26,0) (1) 


Ww =>4 5 
A(9,p)R ~ 26,9) 


Where R array correlation matrix 
R = x(t) * x(t) ) 


and dg is the signal steer vector on array elements can be expressed as [13] 


ag) = Jexp (Gene sin (6) | form =1,2,...,M (3) 
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Thus, the beamformer weights are selected by maximum value, to select the element and minimum 
mean value of output power depending on the number of user inside the coverage while maintaining unity 
response in the look direction [17]. 

The output signal power of the array using optimum weight vector from MVDR beamforming 
method is given by MVDR spatial spectrum for angle of arrival estimated by detecting the peaks in this 
angular spectrum as [18, 19] 


1 
Puvor(@,?) = Gg pe-tae a) (4) 


The standard Capon beamformer has better resolution than the conventional beamformer provided 
that the array steering vector corresponding to the signal of interest is accurately known [20]. 


2.2. LMS 

Least Mean Square (LMS) algorithm depends on continuous update to its parameters. LMS 
algorithm calculates the weight vector using repeated procedure. It makes use of the available data to 
estimate gradient vector. Weights should be in the negative direction of the gradient vector which produces 
minimum square error [21]. The following equations describes the LMS operation [22] 


y(n) = w"(n)x(n) (5) 
e(n) = d(n) — y(n) (6) 
w(n + 1) = w(n) + A. x(n). e(n) (7) 


Where, y(n)=array output, x(n)=input data vector, d(n)=reference signal, w(n)=weight vector, e(n) =error 
signal, and A=step size (0 < A< A,,,,), where Aa, is the largest Eigenvalue of correlation matrix R as in (2) 

Any change to the input data characteristics require changing the filter weights to adapt to these 
changes. The target of the filter is that e(n) becomes zero [23, 24]. 


3. SYSTEM MODEL 

The Massive MIMO system model is represented uniform Rectangular array (URA) with M element 
in the x-axis and N element in the y-axis and spacing between elements d. T different sources transmit signal 
S; at the same time, where 1<1<T. As in Figure 1. 





Figure 1. Uniform rectangular array (URA) 
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The received signal at (m, n) th antenna for "1<m<M,1<n<N" , 1s represented by the following: 


2nd 


: ; 2nd . 
Or S(ieihm Gn 0;+(n-1)(F2) Cos disin g +Wann (8) 


After collecting the received signal in M x N matrix, the final received signal from all sources could 
be represented by two equations using the Separated Steering Matrix method (SSM) in [25] 


X,=A,StW and X,=A,,St+tW (9) 


Where, S=[S,, Sinacsel , W is M x N noise matrix, A,, represents the steering matrix with the elevation 
angle 0, While the A,,, contain both azimuth and elevation angles. 





1 1 = 1 
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The array factor (AF) of M elements in x-dimension can be expressed as, 
AF(@) = "rt Ay Wipe (12) 
And (AF) of N elements in y-dimension can be expressed as, 
AF(g) = yo Aupy, Wuvy (13) 
Normalized Array factor 1s, 
AF 
AP orn ~~ AB ie (14) 


4. PROPOSED METHOD 

The proposed beamforming method combines two adaptive beamformers. This beamformer work in 
cascade form of two stages, where the first stage is the MVDR Beamformer and the second stage is LMS 
beamformer. These two beamformers takes feedback information from the proposed SSM DOA for 
continuous determining and acquisition of users connected to BS. Figure (2) shows a block diagram of 
beamformer. MVDR belongs to the blind beamformers that doesn’t require knowledge of a reference signal 
and, therefore, there 1s no need to adjust the weights of the array. It has one limitation that it requires a prior 
knowledge about the direction of the interference signals. This causes a tangle in apply since the directions of 
interference signals vary over time. 

To overcome this shortage, here comes the role of the second beamformer. The LMS beamformer is an 
adaptive beamformer belongs to the non-blind beamformers and it adaptively varies the antenna array pattern 
to produce nulls toward the direction of interference. 
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Figure 2. Block diagram of hybrid MVDR-LMS 





Hybrid MVDR-LMS beamformer, will benefit from constant initial weights of that MVDR in 
equation (4) to provide fast start in catching the initial user location, while LMS will use equation (7) to 
provide the ability to keep track to user direction and produce null in the interference direction. i.e. this 
proposed algorithm will get use from the good specifications of both MVDR and LMS. The algorithm is 
described by mean of flowchart as shown in Figure 3. 
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Figure 3. Flowchart of hybrid MVDR-LMS 


5. RESULTS AND DISCUSSION 

To show the performance of the proposed hybrid MVDR-LMS algorithm, MATLAB 201 6a is used. 
A massive MIMO model with 128 antennas in the x-dimension and y-dimension was used with 0.5A spacing 
between elements, number of snapshots 50, and carrier frequency 26 GHz. 

Test is carried with four signals received at the array elements with azimuth angle= {-10, 10, 20, 
30} and elevation angle= {10, 20, 30, 40} with 50 snapshots. The direction of arrival was estimated using 
2D-MUSIC using separated steering matrix (SSM) [25]. Figure (4) shows beamforming using MVDR, LMS, 
and hybrid MVDR-LMS for elevation and azimuth angles using 50 snapshots. 
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(c) MVDR-LMS beamformer 
Figure 4. Different beamformers using 50 snapshots 
The figure above shows that all beamformers are able to form peaks in the desired direction, but 
hybrid MVDR-LMS produced narrower beams with power higher than that of MVDR and LMS individually. 


Even when the number of snapshots is reduced to 10 snapshots, it can produce thinner beams with 
maximized power as shown in Figure 5. 
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(c) MVDR-LMS 


Figure 5. beamformers with no. snaps =10 


A test for the simulation time required for beamforming in MVDR, LMS and MVDR-LMS with 
different number of snapshots was made. Table | shows the comparison in time between these algorithms. 


Table 1. DOA and Beamforming Time 


Beamformer Number of snapshots Time in sec 
MVDR 0.811 
LMS 50 0.393 
MVDR-LMS 0.430 
MVDR 0.811 
LMS 10 0.391 
MVDR-LMS 0.410 
MVDR 0.810 
LMS 1 0.385 
MVDR-LMS 0.388 


Increasing the number of snapshots will increase the time required to generate beam form due to 
increase in the time required for training in LMS and MVDR-LMS. While MVDR time to generate beam 
toward each user is not affected with increasing the number of snapshots, since it doesn’t require any 
training. Although MVDR-LMS require computing beam weights using same equation used in MVDR, then 
goes through several iteration to keep track of incoming signal. It still can perform good timing compared to 
traditional LMS. As shown in Table | the time required to compute weight in MVDR equal 0.811 seconds. 
This means MVDR-LMS takes 0.003 seconds more than traditional LMS to generate beams with power 
higher than that of LMS. 
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6. CONCLUSION 

With the increasing number of mobile systems and fast change in their direction, LTE-5G requires 
tracking beamformers with maximized power toward each user. Ignoring interference sources that consumes 
power is also another issue. The available adaptive beamforming algorithms used to produce weights of 
forming and steering the beam to each individual user has some bugs that make it incompatible with 5G. The 
MVDR-LMS BEF presents beams for multi-users with maximum power thin beams towards targets. These 
peaks have power less than that of MVDR by 10.2% and with side lobes power below 100dB. MVDR-LMS 
can presents beam even when the number of snapshots is reduced to 10 snapshots. 

Increasing the number of snapshots will increase the time required to generate beamform due to 
increase in the time required for training in LMS and MVDR-LMS. While MVDR time to generate beam 
toward each user is not affected with increasing the number of snapshots, since it doesn’t require any 
training. Although MVDR-LMS require computing beam weights using same equation used in MVDR, then 
goes through several iteration to keep track of incoming signal. It still can perform good with 3 ms less than 
traditional LMS. 
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